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ABSTRACT 

Aims. We address the task of modeling soft X-ray and optical light curves of X-ray novae in the high/soft state. 

Methods. The analytic model of viscous evolution of an externally truncated accretion o--disk is used. Relativistic effects near a Kerr 

black hole and self-irradiation of an accretion disk are taken into account. 

Results. The model is applied to the outbursts of X-ray nova Monocerotis 1975 (A 0620-00) and X-ray nova Muscae 1991 
(GRS 1124-68). Comparison of observational data with the model yields constraints on the angular momentum (the Kerr param- 
eter) of the black holes in A 0620-00 and GRS 1 124-68: 0.3 - 0.6 and < 0.4, and on the viscosity parameter a of the disks: 0.7 - 0.95 
and 0.55 - 0.75. We also conclude that the accretion disks should have an effective geometrical thickness 1.5-2 times greater than 
the theoretical value of the distance between the photometric layers. 

Key words, accretion disks - binaries - novae: individual (Nova Mon 1975; Nova Mus 1991) 



1. Introduction 

X-ray novae are the low-mass X-ray close binary sys- 
tems (LMXBs) with a re lativistic component, a black hole, or 



a neu t ron star (see, e.g ., Tanaka & Shibazakl Il996t lEsin et all 



I2000t ICherepashchukl 120001) . Much attention has been drawn 
to X-ray novae because they provide black hole candidates 
with the most reliable mass determ inations dCherepashchukl 
l2000t IMcClintock & Remillardll2003l) . The secondary compo- 
nent is a red dwarf, which nearly fills its Roche lobe, and 
transfer rates are estimated to b e as low as 10~ 10 - 10" 12 
MQ>/yr (iTanaka & Shibazakill9 96). Once every several decades, 
the accretion rate onto the central object rises (up to 10~ 8 Mo/yr), 
and the system flares up as a Nova at all wavelengths. 

In a significant number of X-ray novae the X-ray flux 
declines in an exponential fashion after the peak of an out- 
burst, with a characteristic tim e of 30 - 40 days (see, e.g., 
IMcClintock & Remillardl2003l) . The optical flux also decays ex- 
ponentially but about two times more slowly. The X-ray emis- 
sion of the majority of X-ray novae during an outburst can be 
interpreted as the sum of a hard power-law component and a soft 
component produced by a multi-temperat ure black body disk 
with a maximum temperature 0.5 - 1 keV dTanaka & ShibazakH 
119961) . As a rule, during the first several tens of days after the 
peak, the soft component dominates the power-law component 
(this stage is called the high/soft state). After 100-150 days the 
spectrum becomes harder as the power-law component prevails. 
An outburst usually lasts for several months. 

Up to now the mechanism of the outbursts in X- 
ray novae has not been completely understood, but it 



most probably involves an instability of the ou ter disk re- 
gions dKing & Ritterl Il998h ICherepashchukl 120001) . This insta- 
bility, related to hydrogen ionization, was initially proposed to 
explai n dwarf novae, in which accretion occurs onto a white 
dwarf dMever & Mever-HofmeisteHl98ll:[Smakll984t lC annizzo 
Il993h . Subsequently, the conce pt was extended to X- ray novae. 
VanParadijs (1996) followed bv lKing & Ritterid 19981) suggested 
that irradiation of the outer parts of the accretion disk by central 
X-rays has a significant impact on the dynamics of the evolu- 
tionary cy_cJeJn2^a^wvae ; _Sj^hjrr^iation is low in the dwarf 
novae dDubus et alj200lHLasotal200l|) . 

To explain the spectral evolution of X-ray novae, a 
number of models have be e n propose d (see review by 
IMcClintock & Remillardl 120031) . lEsin et al.l (Il997l) proposed a 
scenario involving a change of cooling mechanisms with de- 
creasing accretion rate. According to their model, at the out- 
burst peak and through the subsequent 100-150 days, the disk 
is geometrically thin and optically thick and embedded in a 
hot corona. As the accretion rate falls, the disk, begining from 
the center, gradually becomes a dvection-dominated, geometri- 
cally thick, and optically thin dNarayan et al.l IT9981) . Using a 
succession of sta t ionary models with different accretion rates, 
lEsin et all d 19971, 120001) describe the X-ray and optical light 
curves and the spectral evolution of X-ray nova Muscae 1991 
and X-ray nova Monocerotis 1975. To explain the observed X- 
ray-to-optical flux ratio they considered an irradiation of the 
disks and assumed that 80-90% of the incident X-ray flux is 
thermalized in the disks, which are thicker than given by the 
standard theory. 
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Naturally, time-dependent phenomena in disks call for a dy- 
namical model. Whereas a stationary hydrostatic approximation 
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can be kept to describe the vertical structure of the disk, a non- 
stationary model should primarily address its radial structure. 

A theory of time-depend ent accretion q-disks in bina ry sys- 
tems hasbeen worked out bv lLipunova & Shakural (l2000l) (here- 
after iLSOOl) and subsequently used to describe X-ray and optical 
light curves o f X-ray nova Monocerotis 1975 and X-ray nova 
Muscae 1991 (ILipunova & Shakurall2002l hereafter |LS02). The 
model applies to the first several tens of days after the outburst 
peak, when the spectral state is high/soft, the hydrogen is ion- 
ized in the whole disk, and before the secondary maximum. The 
disk is considered to be externally truncated due to tidal inter- 
actions in a binary. Fitting the observational light curves of the 
two X-ray novae provided estimates of the a-parameter and re- 
lationships between the black hole mass and the distance to the 
systems. Actually, a new method of estimating the ff-parameter 
of accretion disks was proposed, and the or-parameter turned out 
to be rather high, 0.2-0.6. 

In the present work, we develop the model of lLS02l further. 
We include the effect of self-irradiation of the disk and take into 
account the general relativistic effects on the radial structure of 
the disk and on light propagation near a Kerr black hole. The last 
effect increases the degree of irradiation remarkably compared 
to th e plain space-time geom etry, even for a non-rotating black 
hole dSuleimanov et al.ll2007h . 

In Sect. [2] theoretical aspects of the model are considered. 
The modeling procedure and general features of the solution are 
described in Sect. [3] Application of the model to the X-ray novae 
A 0620-00 and GRS 1 124-68 follows in Sect. ED 



2. The model 

2. 1 . Time-dependent accretion disk 

The theory of time-dependent accr e tion h as been put forward in 
the work of lLvubarskii & Shakura (Il987l) and developed for an 
cr-disk in a close binary system by IlsooI 

A diffusion-type differential equation of non-stationary ac- 
cretion is solved analytically assuming specific outer boundary 
conditions. The size of the disk is fixed due to tidal effects, and 
the accretion rate through the outer boundary is negligible com- 
pared with the rate through the inner boundary of the disk. The 
opacity coefficient is a power-law function of density p and tem- 
perature T. Separate regimes with scattering and Kramers ab- 
sorption are considered. The vertical struct ure is treated anayti- 
cally o n the basis of a method proposed by Ketsaris & Shakural 
dl998l) . As a result, the accretion rate at each radius is a power- 
law function of time and of some other parameters. In the case 
of Kramers opacity, it is given bv (iLSOOh : 



M(R,t) = 2n 



1.224^) 



ho 



14/5 



D(t + 5t) 



10/3 



(1) 



Here ho — yjGM\±R oa i is the specific angular momentum of the 
disk matter at the outer radius; 6t has dimensions of seconds and 
normalizes the initial accretion rate; D is a dimensional "dif- 
fusion constant" depending on the <r-parameter, the black hole 
mass, and the opacity coefficient Kq; function y{%) is a polyno- 
mial truncated to three terms: 



y(0 = 1-43- 1.61 £ 5/2 + 0.18<f 



(2) 



where £ = y/R/R ut is the ratio of the angular momentum at 
a given disk radius to the angular momentum at the outer disk 
radius. 



One can infer from (HJ that the accretion rate near the cen- 
tral object hardly depends on R. The disk dynamics (variations 
in the accretion rate, in particular) are defined by its structure 
at outer radii, because it is there that the characteristic viscous 
timescale is the longest. Changes in the accretion rate, occurring 
on a smaller time scale at inner radii, respond to the course of 
evolution at the outer disk. Thus, to analyze the light curves of a 
fully ionized accretion disk in any spectral band, one can use so- 
lution (Q~|i written for the case of Krame rs opacity k - KppT ^ 1 ^ 2 . 

We take k = 5 x 10 24 cm 5 K 7/2 g~ 2 dFrank et alj|2002l chap. 
5), which is the best fit to the real opacity of matter with so- 
lar abundances when hydrogen is fully ionized and bound-free 
absorption by heavy elements is present. Then the diffusion con- 
stant D takes the form 



D = 7.56x 10 3 V /5 Ou/0.5) 



-3/4 



Mbh 



(3) 



The numerical constant in (O incorporates the opacity coeffi- 
cient /to and a factor about 1 that accounts for t he vertical struc- 
ture solution (for more details see lLSO0ULSO2h . 

The key parameters for computing the emission of the disk 
are the effective temperature and the h alf-fhi ckness of the disk, 
crucial for calculating self-irradiation. ILSOOI provides the half- 
thickness of the disk for the zones with Kramers opacity: 



zo(R) 
R 



-1/4 



R 

Rout 



1/20 



3/20 



Here the numerical constant contains a factor about 2.6 that ac- 
counts for the v ertical structure solution and thus depends on 
the value ko (see ISuleimanov et al.l (120071) and Eq. (36) therein). 
Function /(£) is such that y{%) = /"'(£) and 



/(£) = 1.43 £ - 0.46 f' 2 + 0.03 f . 



(5) 



Figure [T] shows the ratio zo/R versus radius obtained by solving 
the vertical structure equations with the tabulate d opacity val- 
ues calculated for the solar chemical abundance (iKuruczll 19941: 
ISuleimanovIl 99 ll ILipunova & Suleimanovl2004l) . It can be seen 
that the half-thickness of the disk calculated from the real (tabu- 
lated) opacity is reproduced satisfactorily by Eq. (@). 

The effective temperature of the disk may be written as 



(6) 



where cr is the Stephan-Boltzmann constant; wk = ^GM/R 3 ; 
<f>(R) is a "structure function" such that for a stationary, not- 
truncated disk in the Newtonian limit <S)(R) — 1 - ^Ri B /R. Taking 
genera l rel ativistic effects into acco unt following lPage & Thornel 
(119741) and lRiffert & Heroldl (1 19951) with zero torque at the inner 
disk radius, we have 



OCR) 



xq 



1.5aln^ 

*0 



A-B-C 



Ax 3 - 3x + 2a) 



Rout m 

R y(0) 



(7) 
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Fig. 1. Relative half-thickness of the disk calculated for Kramers 
opacity (solid line) and for the tabulated values of opac- 
ity (dashed line). The pertubations in the dashed line are caused 
by opacity variations. 
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/ arccos a + n 



)■ 



/arccos a 



xo - 



and Rq = GM\,h/c 2 . The term in the first square brackets 
of Eq. (0 represents a general relativistic correction and ap- 
proaches 1 at large radii. The term in the second square brackets 
of Eq. (O results from the non-stationary solution, and it equals 
approximately 1/1.43 near the outer radius and a; 1 at small 
radii. 

At the outer edge of the disk there are additional sources of 
heating. Some type of collision between the disk and the incom- 
ing matter may cause a development o f a hot area. Moreover, the 
tidal interactions rise steeply at R out dlchikawa & Qsakil 1 1 9941) 
and lead to a corresponding energy release. We believe that nei- 
ther of the above effects contributes significantly to the optical 
flux during the outburst and do not include them in the present 
study. 

2.2. X-ray spectrum 

In the vicinity of a compact relativistic object, the outgoing X- 
ray spectrum is disturbed by the effects of Doppler boosting due 
to the di sk rotation, gravitat ional focusing, and the gravitational 
redshift dCunninghaml 1 9751) . To calculate the spectrum of a disk, 
as detected by an o bserver at infinity, we use a computer code by 
ISpeifh et~ai1 (1 19951) . which provides a transfer function allowing 
for these effects. 

The transfer function is applied to the local X-ray spec- 
tra, which are diluted Planck spectra with a color temperature 
T C (R) = fcTyisiR) due to Compton scattering, whe re / c = 
1.7 + 0.2 (IShimura & Takaharal[T99llDavis et al.ll2005l) . 
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Fig. 2. Upper panel: X-ray spectra of the accretion disk around a 
Kerr black hole calculated for different disk inclinations. Lower 
panel: incident-to-intrinsic flux ratio vs. radius. For both panels 
the Kerr black hole has a = 0.998. The solid lines show the 
impact of the relativistic effects on light propagation and can be 
compared with results obtained in the Newtonian limit shown by 
the dashed lines. 



Figure[2]illustrates the general relativistic effects on the spec- 
trum of the accretion disk around a Kerr black hole (the Kerr 
parameter a = 0.998). Shown are the spectrum as seen by a re- 
mote observer for a disk inclination = 43° (cos# = 0.73) and 
the spectrum for = 87° (cos = 0.05). Naturally, the spectral 
shape of the X-ray emission received by the outer disk (0 is large) 
differs significantly from what is seen by a remote observer. The 
intensity of the illuminating flux is significantly higher than cal- 
culated without general relativistic effects. 

There are different methods of computing general relativistic 
effects on light propagat i on. An example of another implemen- 
tation is one by iLi et all (I2005I) . which is realized in the model 
"kerrbb" in the xspec package. However, kerrbb is restricted to 
inclination angles < 85°. Thus we cannot readily use it to obtain 
the flux illuminating the outer disk, the calculation of which we 
describe in the following section. 

To calculate the X-r ay flux at the Earth we take int erstellar 
absorption into account (Mor rison & McCammonfl983l) . 



4 



Suleimanov, Lipunova & Shakura: Modeling of non-stationary accretion disks in X-ray novae 



CD 



2.0 



1.6 



1.2 



0.8 - 



0.4 - 



0.0 



1 1 1 1 1 1 

■ T =8*10 6 K 




max 




- 1.5-6 keV 




- 


' a=0 




^0.9 




0.9981^ 


i 


i 



0.0 0.2 



0.4 0.6 
COS 



0.8 1.0 



Fig. 3. Function for the different values of the Kerr param- 
eter. See also Fig. 9 in lSuleimanov et all (|2007|) . 



2.3. Self-irradiation of the disk 

The flux falling onto a unit area that is normal to the emission 
outgoing at an angle 6 to the disk axis, can be written as 



\nr L 



(9) 



where Lx is the luminosity of the central X-ray source and r 
the distance of the unit area from the center. In plain space-time 
geometry, for a flat disk as the central source, ¥(0) = 2 cos(fl). 
We us e the function ¥(0) provided by the code of ISpeith et al.l 
(119951) . Examples of "¥(0) for some values of the Kerr parameter 
are shown in Fig. [3] 

To calculate flux absorbed by a unit area of the disk, 
one should all ow for not all of the incident flux be- 
ing thermalized dShakura & Sunyaevlll973l : Ivan Parachjsl 1 19961 : 
lKing&Ritterlll998l) : 



F m (R) = rj 



4-ttR 2 



dz Zq 
dR R 



(10) 



where 77 is a thermalization coefficient, and the term in square 
brackets stands for the inclination of the disk surface to the in- 
cident rays. Assuming the disk to be geometrically thin, the dis- 
tance r is replaced by the radius of the disk R at which the area 
in question is located. 

The geometrical factor in the square brackets of Eq. ( fTOt is 
remarkably different for the cases of the time-dependent model 
and a stationary disk. A stationary disk with M = M(0, f) has ap- 
proximately the same zo/R as in Eq. ©. However, value dzo/dR 



differs significantly. Using Eq. © we obtain the geometrical fac- 
tor for a non-stationary, externally-truncated disk, 



dzo 
dR 



ill 
R 



ill 
R 



20 40 V *out /(£), 



(11) 



which is « 1/20 (zo/R) for R « R out , while the geometrical fac- 
tor for a steady-state disk model (Shakura & Sunyaev 1973) is 
1 /8 {zo/R)\ thus, the time-dependent disk near the outer radius is 
less irradiated. 

The factor 77 can be thought of as a coefficient for the re- 
processing of incident X-rays ("reprocessing efficiency"), i.e. 
the fraction of X-rays absorbed in the under-photosphere lay- 
ers. Values of 77 suggested by researchers differ significantly. For 
low-mass X-ray binaries, Ide Jong et al.l (1 19961) obtained a low 
value for reprocessing efficiency 77 ~ 0.1. In contr ast, mode- 
ling th e outbursts of GRS 1124-68 and A 0620- 



-OO. lEsin et al.l 

(120001) assume that almost all the incident X-ray flux is thermal- 
ized, 77 = 0.8 - 0.9. 

Taking another approach, a value for the reprocessing ef- 
ficiency can be derived from an accurate consideration of the 
atmos phere of an accretion d isk exposed to the incident radi- 
ation ( Suleimanov et a 1. 1999). For a disk with ionized hydro- 
gen, a decisive factor is the spectrum of the incident X-rays. 
Only rather hard X-rays, of energy E > 2 keV, penetrate the 
layers where the optical continuum is produced. Softer X-rays 
are abs orbed higher in the atmosphere and reprocessed to the 
far U V dSuleimanov et alJll999tlFT-KJiourv & Wickramasinghel 
119991) . To raise the optical emission from the disk, the tempera- 
ture should increase at the depths where the optical continuum is 
formed. However, most of the X-rays with energies over 10 keV 
are reflected due to electron scattering. Thus, it is mainly X-rays 
in the interval 2- 10 keV that are reprocessed to optical emission. 
The high/soft spectral state of X-ray novae is characterized by 
soft X-ray spectra with a correspondin g black body temperature 
~ 0.5 - 1.0 keV. Using Eq. (22) from lSuleimanov et alj d 19991) 
we obtain low v alues o f the reprocessing efficiency, 0.1 - 0.2. 

lEsin et all (120001) generally conclude that the outer 
parts of accretion disks in the high/soft state should be 
si gnificantly thicker than fo llows from the standard model 
of IShakura & Sunvaevl d 1973b . Concerning this result we bring 
up the following argument. The irradiation cannot increase the 
thickness of the disk beyond the values predicted by Eq. ©, 
either in the regions where Kramers opa city applies, or in 
the regions with partly ionized hydrogen dDubus et al.l 119991: 
ISuleimanov et al.l 120071) . However, a layer of hot coronal gas 
could be present abo ve the disk due to the irradiation of the outer 
disk b y hard X-ravs (IBegelman et ail 1 9831: iBegelman & McKe" " 



1981" 



Ko & Kallmanl 1 1 994 iRozanska & Czernvl \199t 
Proga & Kallmanl 120021: Uimenez-Garate et al.l |2002|) . Such 
a corona, optically thick in the radial direction, effectively 
increases the geometrical thickness of the disk and augments 
the degree of self-irradiation. Furthermore, the corona could 
be inhomogeneuos and contain optically th ick gas blobs as 

[Toposed for supersoft X-ray sources by ISuleimanov et al.l 
20031) . In this case the reprocessing efficiency can increase 
significantly due to multiple scattering between the blobs. 

There is another possibility for increasing a disk's thick- 
ness. Recent simulat i ons of magnetorotational turbu lence (e.g. 
iMiller & Stond2000tlTurnerl2004lHirose et al.l2006l) show that 
magnetic pressure is significant in the upper layers of the disk, 
probably resulting in a vertically extended atmosphere. Thus, the 
optical-depth unity surface can be located a fac tor of two above 
the "standard" disk surface (Hirose et al . i 2006l) . 
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In view of the present uncertainty, we consider the repro- 
cessing efficiency as a parameter. For the same reason, we do 
not calculate the half-thickness of the disk using the tabulated 
values of opacity coefficient but utilize Kramers approximation. 
In Eq. ( fTOt . the presence of the corona is treated as the increased 
thickness of the disk, which equals zo multiplied by a chosen 
factor. 

It should be noted that consideration of the hard spectral 
component, which is emitted by a central source (possible cen- 
tral corona) and characterized by ^{9) ~ 1, would enhance the 
self-irradiation of the disk. 



2.4. Optical emission from the disk 

The spectrum of the time-dependent disk at any moment of time 
is calculated by integrating local spectra over radius: 



r Ro, 



B A (T)RdR. 



(12) 



At the outer parts of the disk the local spectra are Planck func- 
tions with an effective temperature T increased due to the inter- 
cepted X-ray flux 



(13) 



where r v ; s is the effective temperature (Eq. |6| in the absence of 
X-ray irradiation. The spectral flux density at the Earth is given 
by 



h = T~F cos ' ex P( _T -i) > 
2nd- 



(14) 



where d the distance to the system, i the inclination of the disk, 
Ta is the optical depth for interstellar absorption. 

The unabsorbed optical and UV spectrum of a self-irradiated 
disk is shown in Fig. [4] 



3. Modeling procedure 

The physical parameters of X-ray novae (mass, Kerr parameter, 
viscosity parameter a, etc.) are sought by comparing modeled 
and observed light curves. A particular model of the disk is de- 
fined by a set of input physical parameters. Some of these can be 
varied and some are assumed to be known from observations. 

To compare the model with observations, we chose several 
values calculated in a specific way from observational data- "fit- 
ted light curve parameters"- and fit them by the model. 



3.1. Fitted light-curve parameters 



As shown in ILS021 the X-ray and optical light curves of the X- 
ray novae considered are represented well by linear functions in 
"log(flux)-time" coordinates throughout the first 40-60 days af- 
ter the peak. Thus each light curve can be described by a pair of 
parameters: the slope and the flux at a reference time, which are 
to be fitted. This approach allows one to significantly reduce the 
computer time required for finding satisfactory model parame- 
ters. 

In the course of the modeling, we reproduce the following 
values: the slope of the X-ray light curve ax, the slope of the 
optical light curve in the B band as, the X-ray flux Fx, the flux 
Fb, the color index B - V, and, optionally, U-B. All parameters, 
except for the first two, are calculated for a reference time. 
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Fig. 4. Optical and UV spectrum of a modeled disk of 
A 0620-00 with self-irradiation (the solid line) and without (the 
dashed line). Vertical lines denote the positions of the centers of 
B and V bands. 

3.2. Input parameters of the model 

The free input parameters are: the black hole mass Mbh, its Kerr 
parameter a, the viscosity parameter a of the disk, the time pa- 
rameter 6t to fix the accretion rate at time t — 0, and the re- 
processing efficiency rj. Other variable parameters of the model 
are the color-to-effective temperature ratio / c and the ratio z/zo, 
where z is the height of X-rays' interception. 

The parameters that are more or less fixed include: the op- 
tical component mass M , the mass function f(M), the period 
of the binary P, and the interstellar absorption, expressed by pa- 
rameter Nh for an X-ray band (the number of H atoms along the 
line of sight per cm 2 ) or by the color excess E(B - V) for optical 
bands. 

The mass of the optical component M and the interstellar 
absorption are not very precisely known. In lLS02l thev were var- 
ied within narrow limits, and it was shown that the effect of the 
variation was small, so we fix their values. 

Given the binary mass, one obtains the semi major axis of 
the binary system 



a = 



G (M b h + M ) P* 
4n 2 



2\l/3 



(15) 



and can calculate the outer radius of the disk approximately. We 
take this equal to 80% of the radius of a sphere whose vo lume 
equals that of the black hole's Roche lobe (Eggl etonll983l) : 



0.8 5 



0.49 q 



2/3 



0.6 q 2 ' 3 +ln(l + 9 1/3 ) 



(16) 
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where q - Mhh/Mn- Thi s approximation agrees with calcula- 
tions by iPaczvnskil ( I 19771) and with the direct estima tes of disk 
sizes in cataclysmic variables dSulkanen et al.|[l98ll) . One can 
also derive the inclination i of the binary plane: 



arcsin 



( f(M)(M hh + M ) 
Ml 



2 \ 



1/3 



(17) 



It is necessary to mention that the inner disk can have another in- 
clination to the line of sight. The black hole angular momentum 
can be misaligned with that of a binary due to asymmetric previ- 
ous supernova explosion. In this case the inner disk will tilt to- 
ward to equatorial plane of the black hole ( Ba rdeen & Pettersonl 
1 1975b . We do not consider this possibility in the present study 
and assume that the inner and the outer parts of the disk have the 
same inclination; therefore, some of our results can be biased, 
black hole angular momentum in particular. 

The distance to the system is defined in the course of the 
modeling. As the limits on black hole masses and distances are 
known from observations, they are used to restrict the allowed 
ranges of other parameters. 

3.3. Qualitative investigation of nonstationary accretion in 
a-disks 

Before proceeding to the model results, we discuss the qualita- 
tive features of the nonstationary model. 

An estimate of R out can be derived from Eqs. ( fTBI ) and dT6b 
in the approximation that M D h » M 



flout <* M l b l 3 pV3 



(18) 



Using this relation and Eq. ([TJi, we obtain 



M oc 



M; 



20/9 „25/9 
bh 



a*>Ht + St) 10 ' 3 ' 
Next, the bolometric luminosity of the disk can be written as 



(19) 



L bo i = s(a) Mc 2 7r cr | — | R x . 



(20) 



where s(a) is the efficiency of conversion of gravitational energy 
to radiation, dependent on the dimensionless Kerr parameter a; 
Tx is the maximum color temperature at radius 



Rx~R ln K{3(a)M hh . 



(21) 



Here B{a) represents a dependence on the Kerr parameter. For 
a Schwarzschild black hole a = 0, e — 0.057, B - 6 and for 
a maximally rotating Kerr black hole a = 1,6 = 0.42, B 
I . Recall that the color temperature is / c times greater than the 
effective temperature. From Eqs. (fT9l and (|20b we have 



„l/4 A/'/ 18 P25/36 
rp _ bh 

X °V/ 2 a^it + St) 5 ! 6 



(22) 



lEbisawa et al. (1994) obtained the maximum color temperature 
variations for GRS 1124-68. Figure [5] shows good agreement 
between the observed decrease and the theoretical dependence 
on time (Eq.l22l. 

The fading of an X-ray light curve is faster than that of a 
bolometric light curve, expressed by Eq. ( fT9l . An interval of the 
spectrum that we consider can be approximated by Wien's law 
with a temperature close to the maximum color temperature of 
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Fig. 5. Observ ed maximum color temperature of the disk in 
GRS 1 124-68 (Eb isawa et al.ll 99*41) vs. time in comparison with 
Eq. d22l) where 5t = 105 days expressed in seconds. 



the disk. Thus, variation in the X-ray flux of a soft X-ray tran- 
sient is governed by the decreasing color temperature and the 
associated shift of the Wien-like spectrum to lower energies; a 
well-kno wn exp onential dependence of the light curve on time 
emerges (iLSOOt) . The optical light curves have flatter slopes be- 
cause the main decrease in the disk luminosity happens at high 
energies. 

Let us calculate the slope of an optical light curve, for ex- 
ample, in the B band. Because the optical flux from the disk is 
mainly determined by the reprocessed X-ray emission from the 
central source, we can write 



(23) 



where index k > 1 comes from a relation between the optical 
flux and the effective temperature. For instance, in the Rayleigh- 
Jeans part of the spectrum k — 1 . The time dependence of the ef- 
fective temperature at the outer edge of the disk can be obtained 
from Eq. (flDl : 



^bol / Zo 

R 

or, using Eqs. ©, (Tj9]i, and d20l l 

r 11T <x( £7 ,) 1/4 



oT ilT oc 77 — ^— — — 
R„,„ \ flout 



M w* P uns 

bh 



a ll/12 (f + ( j f) 13/12 

The slope of the B light curve at time t is 
dig L B (t) 13£lge 



a B 



dt 



12 t + 6t 



(24) 



(25) 



(26) 
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Table 1. Observed input parameters and fitted light curve 
parameters". 



Parameter 


A 0620-00 


GRS 1124-68 


P [days] 


0.322 1 - 2 


0.433 3 


KM) 


2.7M^ 2 


3.0M^ 


M a 


O.5M 


O.SMq 


N H [atoms cm" 2 ] 


3 x 10 21 b - c 


1.4 x 10 21 c - 6 


E(B - V) 


0.35 3 


0.3 C 


ax 


-0.0149 ± 0.0004 


-0.0125 ± 0.0005 


F x 


100 photons s -1 cm -2 


6.8 x 10" 8 erg s~' cirr 2 


a B 


-0.0074 ± 0.0003 


-0.0060 ± 0.0004 


F B [erg s cm 2 ] 


(1.24 ± 0.07) x 10- 10 


(1.12 ± 0.03) x 10-" 


B-V 


0^24 + (F03 4 


(T27 ± OW 4 


U-B 


-OT75 ± 0?05 





Notations are explained in Sects. [3TTl and l3,2l Absolute fluxes are 
derived for the 30th day after the outburst peak. The slopes of 
the X-ray light curves and the X-ray fluxes correspond to differ- 
ent X-ray intervals: 1.5 - 6 keV for A0620-00 and 1.2-6 keV for 
GRS 1124-68. 

Interm edia te between the v alues obtained by IWeaver & Williamsl 

( fl974l) and lWu et al.1 ( fl983l) . 

N H is obtained fr om E(B - V) fol l owing IZombeckl JT990b. 



Refe rences. (1 ) 
(3) IWuet al 



lOrosz et al 



CherepashchuU feOOd); (2) IChen et all {19971) : 
.19831) ; ( 4) iLipunova & Shakurd 112003); (5) 
1996»: (6) Idella Valle et all dl99lh : ICheng et all 



dl992h ; IShrade7& Gonzal ez-Riestral ( 1 19931) 



and is essentially determined only by the parameter 6t and vice 
versa. Our numerical estimations yield k w 1.9 for the spec- 
trum of an irradiated disk. Substituting values of slopes ag from 
Table [TJ into Eq. (26) , we obtain that 5t is approximately 90 and 
100-110 days for A 0620-00 and GRS 1124-68 (for the 30th 
and the 50th days after the peak, respectively). 

The color index B-V obviously determines the temperature 
at the outer disk edge T m (R = R out ). Thus, the parameters a, M\,h, 
and a, which produce a satisfactory model, should be related to 
meet the condition T m (R = R ou{ ) ~ const, defined by Eq. ( 1251 ). 
To conclude, we note a basic feature: for a constant black hole 
mass and 7x, an increase in the Kerr parameter will cause the 
distance to the binary to decrease (cf. Eqs. d20b and d2Tl)). 



4. Outbursts of Nova Mon 1975 (A 0620-00) and 
Nova Mus 1991 (GRS 11 24-68) 

Observational parameters for Nova Mon 1975 and Nova Mus 
1991 are summarized in Tab le [D We take B and V opti- 
cal light curves collect ed in lLS02l The U light curve o f 
A 06 20-00 is taken from iLiutvil (1 19761) . Ivan den Berghl dl976l) . 
and iDuerbeck & Walter! ( 19761) and trans ferred to units of 
erg cirr 2 s in the same way as described in ILS02L 

By and large, the outbursts are very much alike. Let us point 
out the distinctions. In GRS 1124-68, all light curves decline 
more slowly: the decay time scales, or "e-folding times", of 
the X-ray and optical light curves are 35 days and 68 days for 
GRS 1 124-68 as against 30 days and 57 days for A 0620-00. In 
addition, the optical-to-X-ray flux ratio at the peak in A 0620-00 
is ~ 2 times higher than in GRS 1124-68. The color index 
(B-V) is greater for GRS 1 124-68 (the object is "redder"), al- 
though the color excess E(B - V) is higher for A 0620-00. 

The input parameters of the model are varied within the fol- 
lowing limits: M bh : 5-15M©, the Kerr parameter a: 0-0.9981, 
the ff-parameter: 0.1 - 1, the time normalization parameter 6t: 
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0.1 



5t = 90' 
d = 1 .2 kpc 
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Fig. 6. Comparison of the ob served spectrum of 
A 0620-00 ( Long & Kestenbaum|[l978l) and the modeled one 
for t = 50 d . The solid line represents the case when interstellar 
absorption is taken into account with A^h = 3 x 10 22 atoms cm -1 ; 
the dashed line shows the unabsorbed modeled spectrum. The 
lower panel shows a enlarged area from the upper panel. 



80 - 140 days. The modeling results critically depend on the 
reprocessing efficiency 77, which we vary between and 1 . 

4.1. A 0620-00. Importance of irradiation 

We use the 1.5 — 6 keV light curve obtained with the SAS-3 
CSL-A dKuulkersI [19981) . The flux data are converted from 
counts s~' to photons cirr 2 s~' assuming a Crab-like spec- 
trum. The fo llowing normaliz ations are adopted: 1 Crab * 
136 counts/s dBuff et al.Hl977[) . 1 Crab * 5 photons cm" 2 s 1 
in 1.5-6 keV. We believe the resulting error to be acceptable 
because 1.5-6 ke V and 3-6 keV fluxe s in Crab units are 
close to each other (K aluzienski et al.ll 1977b . and thus the spec- 
tral shapes of A 0620-00 and the Crab should be similar in the 
X-ray range in question. In Fig. [6] we compare the X-ray spec- 
trum calculated from the m odel with the spectrum observed by 
lLong & Kestenbauml(ll978l) . The latter was obtained on the 70th 
day after the peak, when the flux had the same magnitude as on 
the 50th day dKuulk ers 1998). We conclude that the chosen nor- 
malizations are satisfactory. 

During the interval of 20 - 40 days after the peak the multi- 
temperature black body disk component dominates the spectrum 
(see lLS02h . In Table [T]the slopes of the light curves are calcu- 
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0.0 0.2 0.4 0.6 0.8 

Kerr parameter 



Fig. 7. Resulting parameters of A 0620-00 vs. Kerr parameter. 
Models with f c = 1.7 are denoted by circles (z = zo and rj-l), 
stars (z = 2zo and 77=0. 5), and pluses (z = 2zo and 77=0.5 with 
limb darkening). Models with / c = 1 are denoted by diamonds 
(z = Zo and 77=1) and triangles (z = 2zo and 77=0.5). Models 
den oted by filled stars ha ve masses and distances in agreement 
with lGehno eTail d2001al) . 



lated for this time interval, and values Fx and Fb correspond to 
the 3 0fh day . 

In lLS02l the difficulty of reproducing the slope of the optical 
light curves of A 0620-00 was revealed. Including reprocessed 
emission helps to solve the issue. In the case of the self-irradiated 
disk, the slope of the optical light curve steepens, partly because 
the incident flux declines faster with time than the intrinsic flux, 
but mainly because St has to be decreased to adjust the optical- 
to-X-ray flux ratio, and 5t is uniquely related to the slope (see 
Eq. d26|i). Modeling results are presented in Fig. [7] The relation 
between allowable values of the black hole mass, the Kerr pa- 
rameter, and the a-parameter reflects the constancy of the value 
Tx- An example of the modeled and observed light curves is 
shown in Fig. [8] 

There are two parameters that define the effect of irradia- 
tion (see Eq. (TlOb): the reprocessing efficiency 77 and the ratio 
z/zo- In the Newtonian limit the combination 77 (z/zo) 2 is deci- 
sive, but in general relativity this is correct only approximately 
due to the relativistic function ^¥(9). Our calculations show that 
the observed parameters can be fitted with r](z/zo) 2 ~ 1 - 2. 
Evidently, the reprocessing efficiency 77 can hardly be as high as 
1. Therefore, the accretion disk in this X-ray nova must have 
an effective geometrical thickness greater than that given by 
Eq. ©. One can introduce the so-called irradiation parameter 



C t o rewrite Eq (TTOTt : F in (R) = C(L x /4jtR 2 ) (see, for exam- 
ple, iDubuTeTayllOOll). For the model presented in Fig. [8] in the 
time span 10-50 days, the irradiation parameter, calculated at 
the outer radius R out = 1.59 x 10 11 cm, changes in the range 
~ (8 - 6) x 10~ 4 . 

There are no satisfactory models with the Kerr parameter a > 
0.7 for / c = 1.7 (Pig. O. We also investigate the influence of 
parameter f c . Results with f c = 1 are presented in Fig. Qby 
diamonds (for 77 = 1 and z/zo = 1) and triangles (77 = 0.5 and 
z/zq = 2). One can see that models with f c = 1 have higher 
values of the Kerr parameter, which compensate the decrease in 
the color temperature. Other parameters (the black hole masses, 
ff-parameter, and the distance), though somewhat decreased, are 
still in the same range. Therefore, any uncertainty in f c leads 
mainly to a deviation of the Kerr parameter. 

In some of the models, the effect of limb darkening was in- 
cluded. For the local X-ray spectra, the limb darkening was taken 
into account assuming the dominant electron scattering, and the 
grey atmosphe re (Eddington) approach was used for the local 
optical spectra dChandrasekhadl 1 9501: ISobolevll 1 949l 1 1 9691) . The 
corresponding results are shown in Fig. Q by pluses. Taking ac- 
count of limb darkening leads to some increase in the black hole 
mass, the a-parameter, and the distance to the binary, not altering 
the results qualitatively. 

Independent estimates of the system parameters can be used 
to constrain our results further. Available in th e literature are 
estimates of the distan ce d = 1+0.4 kpc dShahbaz et al.l 
ll994tlBarretetalJll996l) a nd the black hole mas s: M bh = 10 ± 
7Mg, dShahbaz et al.ll 1994 . M bh = 4.6 ± 1.1M© dHaswell et al.l 
1993 ). The h i ghest q uoted accuracy is of the estimates given by 
iGelino et al.l feoOlab : M bh = 11 ± 1.9M© and d = 1.116 ± 
0.114 kpc. The constructed models agree with these values if 
a = 0.7 - 0.95 and the Kerr parameter a = 0.3 - 0.6 (Fig. [7j. 

4.2. GRS 1124-68 



lEbisawa et al.l d 19941) fit the spectra obtained during the 1991 
outburst of GRS 1 124-68 with a model comprising two compo- 
nents: a hard power-law c omponent and a soft component emit- 
ted by a multicolor disk dMitsuda et alJI 1 984b . They present a 
time evolution of model parameters: the maximum temperature 
and the inner radius of the multicolor disk, as well as the photon 
index and the hard 2-20 keV flux from the power-law com- 
ponent. This information allows us to calculate the unabsorbed 
1.2-6 keV flux. In Table [TJ the slopes of the light curves are 
calculated for a 20-50 day interval. 

It is possible to fit the light curves with plausible values of 
reprocessing efficiency, 77 < 0.5, keeping the standard thickness 
of the disk. For this system the relation 77 (z/zo) 2 ~ 0.5 - 0.6 is 
fulfilled. One gets almost same results with three different sets 
of parameters: (1) z = Zo, tj = 0.5; (2) z = l.5zo, H = 0.25; 
(3) z = 2zo, 77 = 0.15 (see Figs. |9l and ITOb. The parameter 5t 
lies in the range from 100 to 115 days in all models. Figure [9] 
presents allowable intervals for the mass, the Kerr parameter, 
the ^-parameter, and the distance, which lies in the interval from 
2.5 to 6 kpc. 

Figure [TTJ illustrates the modeled light curves for parame- 
ters M bh = 8M©, a = 0.55, a = 0, 8t = 115 d , 77 = 0.25, 
and z = 1.5zo- During the interval of 10 - 60 days, the irradi- 
ation parameter C (see Sect. 14.11 ). calculated at the outer radius 
R out = 1.62x 10 11 cm, decreases from ~ 3 x 10~ 4 to ~ 2.5 x 10 -4 . 
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Fig. 8. Comparison of the modeled and observed light curves of 
A 0620-00. The U flux is shifted upward by 0.2 . 

4.2.1. Distance to GRS 1124-68 

Independent estimates of the distance and the bl ack hole mass 
are: d = 3.4 + 0.6 kpc and M hh = 7.5 ± 3.5M© dShahbaz et al 
1997I) . d = 5.5 + 1 kpc and M bh = 5 - 7.5M dOrosz et al 



1996 ). The high e st accur acy f or the black hole mass is claimed 



bv iGelino efatl d2001bl) and iKhruzina et all d2003l) : M bh = 
6.95 + 0.6Mp and M hh = 7.15 + 0.45 M©, respectively. 
IGelino et alj d2001bl) estimate the distance as 5.1 kpc or even 
5.9 kpc (lGelinoll2004l) . 

The c onstructed models ca nnot agree with th e distances ob- 
tained bv lGelino et al.1 d2001bl) andlGelinol (|2004), tending more 
to the estimates of lShahbaz et alJdl997l) . For the black hole mass 
M = 7.5 + 3.5M©, z = (1.5 - 2) Zo , and tj = 0.15 - 0.25, we 
come to the following allowable intervals of parameters: the dis- 
tance 3-4 kpc, the Kerr parameter a < 0.4, and the a-parameter 
0.55 -0.75 (Fig. 

The above estimates of the distance are derived from an anal- 
ysis of light curves and spectra in optical and IR bands ob- 
tained during the quiescent state of the X-ray nova. The ba- 
sic element for determining the distance in quiescence is the 
spectral class of the secon dary. It is found from optical spectra 
(~ K4V, lOrosz et aTll 19961) and yields a corresponding effective 
temperature (~ 4500 K). However, the optical light curves of 
the seconda ry imply a reflection, and it is necessary to assume 
~ L opt dKhruzina et al.l 120031) . Thus, the apparent spectral 
class of the secondary is biased; if the actual effective tempera- 
ture of the secondary is lower, then the binary is closer. 

Another way to constrain the distance is to study X- ray spec- 
tra obtained during the outburst. lEbisawa et alj d 19941) obtain a 
non-relativistic relation between the black hole mass and the dis- 
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All models have f c = 1.7. 



tance, which can be rewritten as follows: 



1/2 



kpc. (27) 



This relation appears to agree with our results (Fig. |5), keeping 
in mind that slight changes due to relativistic corrections might 
be present. 

It is worth noting here the possibility of increasing a dis- 
tance obtained from X-ray data. The inner disk radius can be 
greater than the radius of the last marginally stable orbit ,R ms . 
This ph enomenon is discussed for other Galactic black hole can- 
didates dGilfanov et al.lll999t loiatv et alj|2003l) . The inner re- 
gion, evacuated by the standard disk, becomes hot and optically 
thin, so not contributing to the soft X-ray flux. We tested a num- 
ber of models with the inner radius moved away up to 10/? ms . 
The imposed boundary condition was that the viscous stress ten- 
sor becomes zero at R ms and is non-zero at the inner radius R^. 
We came to the conclusion that all such models, which satisfac- 
torily fitted the data, had higher values of the Kerr parameter and 
unchanged values of the mass and the distance. 

5. Discussion and conclusions 

As shown, the model with non-stationary accretio n onto a com- 
pact binary component, developed for an a- diskin|LSlJ can ac- 
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Fig. 10. Resulting relationships between the mass of the black 
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count for the X-ray and optical light curves of the two X-ray no- 
vae considered, A 0620-00 and GRS 1124-68, in the high/soft 
state. Dedic ated mo deling of the light curves has previously been 
attempted in lLS02l In the present study we extend the modeling 
by including the effects of general relativity on light propagation 
near black holes and self-irradiation of the disks. 

One needs to assume a rather large contribution of repro- 
cessed flux in the optical light curves. A thin disk encounters 
difficulties here in explaining the required amount of intercepted 
flux. A twisted disk could be suggested here as a solution be- 
cause :Jtej^gjojis facing the center under go increased illumina- 
tion dDubus et al.l 1 19991 lEsin et al.ll2000l) . but the total benefit, 
calculated for the whole surface of one disk side, is ambi guous 
due to the effect of self-shielding (ISuleimanov et al.ll2.Q07l) . 

In the current study we prefer another explanation. We sug- 
gest the presence of matter above the disk in the form of a hot 
corona with temperature ~ 2 x 10 6 K ( e.g., corresponding to 
the at mosphere or the "warm corona" of iJimenez-Garate et al.l 
|2002|) . By intercepting X-rays, the corona effectively increases 
the disk thickness. This idea manifests itself in the computation 
of the optical emission when we substitute the half-thickness of 
the disk zo by the value z ~ (1.5 - 2)zo representing the height 
of X-rays' interception. Furthermore, if the corona is inhomoge- 
neous, i.e. consists of cooler blobs surrounded by hotter plasma, 
then multiple scattering between the b lobs could explain the hi gh 
values of the reprocessing efficiency (ISuleimanov et a l. 2003). 
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Fig. 11. Comparison of the modeled and observed light curves 
of GRS 1124-68. 



The outer parts of the disks in A 0620-00 and GRS 1 124-68 
have remarkably different properties. The disk in GRS 1 124-68 
turns out to be thinner and/or and weakly thermalizes the inter- 
cepted X-ray flux. Probably the illumination of the outer disk 
in A 0620-00 is more pronounced because of the higher Kerr 
parameter. If we freeze the disk thickness at a standard value 
and re-calculate models with the parameters used in Figs. [8] 
and QT| we find that the irradiating flux in A 0620-00 is twice 
that in GRS 1124-68. This could be the reason for the more 
developed corona above the outer disk in A 0620-00. Another 
possible source for the extended atmosphere of the disk is the 
magnetic pressure supporting the outer layers of the disk (e.g. 
iHirose et~a l. 2006). This would agree with the higher value of a 
in A 0620-00, given that the magnetorotational instability gen- 
erates the viscosity. 

To bound the resulting parameters we need to invoke in- 
dependent estimates of the mass and the distance. Then for 
A 0620-00 we obtain the intervals: a = 0.7 - 0.95 and the 
Kerr parameter a = 0.3 - 0.6. For GRS 1124-68 we find a = 
0.55 - 0.75 and the Kerr parameter a < 0.4 if the distance to the 
system lies between 3 and 4 kpc. If we do not apply the knowl- 
edge of the mass and the distance, we are nevertheless able to 
confirm the high value of the turbulent parameter: a > 0.4 (mid- 
dle panels in Figs.|7]and|9j. Such values of a are in aproximate 
agreement, but slightly higher than the commonly pr oposed val- 
ues fo r thin, fully ionized disks: a ~ 0.1 - 0.4 (e.g. jKing et al.l 
120071) . 

The accuracy of the estimates of a obtained suffers from 
limitations in the calculations of the X-ray flux made in the 
model and from the observational data. In this regard, we should 
also mention a simplified account of the tidal effects in the an- 



Suleimanov, Lipunova & Shakura: Modeling of non-stationary accretion disks in X-ray novae 



11 



alytic model of a non-stationary disk dLSOOl) . Presumably, the 
accretion disk in an X-ray nova binary has a fixed outer ra- 
dius, at which a delta function-type transport of angular mo- 
mentum to the orbital motion occurs due to tidal interactions. 
A physically more compre hensive consideration, involving a 
real form of the tidal torque dPapaloizou & Pringle 1977tlSmakl 
Il984tllchikawa & Osakilfl 994t lHameurv & Lasotal(200l . could 



revise the solution (Eq. Q]) and thus change the estimates of 
a (see Eq. (T%). From this point of view, the value of a ob- 
tained for the outer disk is a "characteristic" parameter of the 
model. However, an advanced treatment would not change sig- 
nificantly the theoretical M, and thus a, if the tidal interactions 
are small in the whole di sk and rise very steeply only near R out 
dlchikawa & Osakill 19941) . 

Our estimates of the Kerr parameter basically accord with 
a variety of estimates recently obtained for several stellar mass 
black holes: a = 0.77 + 0.05 for M33 X-7, a = 0.65 - 0.75 for 
GROJ1655-40,a = 0.75 - 0.85 for 4U 1543-4 7, a = 0.98 - 1 .0 
for GRS 1915+105, a < 0.26 for LMCX-3 dLiu et alj|2008t 
IShafee et al.l2006tlMcClintock et al.l2006tlDavis et al.l2006h . In 
these works, X-ray spec tra are fitted by spectral models of a rel - 
ativistic accretion disk dLi et alj|2005t iDavis & Hubeny||2006l) . 
and the black holes masses and orbital inclination angles are 
well known. Also, the Fe line method has been used to ob- 
tain an esti mate for the blac k hole spin a = 0.86 ± 0.05 in 
GX 339-4 dMiller et al.ll2008l) . 

The model presented here can only be applied during the 
first few tens of days after the peak and fails to explain the late 
outburst phases. In the low/hard state, the thin disk model can- 
not adequately account for the observed spectra and energetics 
of X-ray nov ae. The ADAF model is believed to be appropri- 
ated instead dEsin et all 1 19971 l2000b . On the other hand, com- 
plex X-ray light curves at the early outburst phases, for example 
those of GRO J 1655-40 and XTEJ1550-564, can hardly be in- 
terpreted by means of the present method. These binaries are 
highly inclined to the line of sight, and the flat thin disks are 
obscured by the outer thick d isks during the early outburst (see 
iNaravan & McClintockll2005l) . The law of accretion rate evolu- 
tion (Eq. [1} also does not hold for the events involving an addi- 
tional enhanced mass transfer from the secondary. 

Our choice of the two X-ray novae stems from their exhibit- 
ing the most classical shape of X-ray nova light curves. The 
prospect of a future determination of a in ionized accretion disks 
in binary systems is likely to come from modeling the evolution 
of disk parameters during the time when the disk spectral com- 
ponent dominates in the spectrum of an exponentially decaying 
X-ray nova in the high/soft state, and preferably with known pa- 
rameters and orbital elements of the binary. 
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